Full orbit simulations of collisional impurity 
transport in spherical tokamak plasmas with 
strongly-sheared electric fields 

C G Wrench^ E Verwichte^ and K G McClements^ 

^ Centre for Fusion, Space and Astrophysics, University of Warwick, Coventry, CV4 
7AL, UK 

^ EURATOM/CCFE Fusion Association, Culham Science Centre, Abingdon, 
Oxfordshire, 0X14 3DB, UK 

E-mail; e.g. wrenchSwarwick .ac.uk 

Abstract. The colhsional dynamics of test impurity ions in spherical tokamak 
plasmas with strongly-sheared radial electric fields is investigated by means of a test 
particle full orbit simulation code. The strength of the shear is such that the standard 
drift ordering can no longer be assumed and a full orbit approach is required. The 
effect of radial electric field shear on neoclassical particle transport is quantified for 
a range of test particle mass and charge numbers and electric field parameters. It 
is shown that the effect of a sheared electric field is to enhance the confinement of 
impurity species above the level observed in the absence of such a field. The effect 
may be explained in terms of a collisional drag force drift, which is proportional to 
particle charge number but independent of particle mass. This drift acts inwards 
for negative radial electric fields and outwards for positive fields, implying strongly 
enhanced confinement of highly ionized impurity ions in the presence of a negative 
radial electric field. 
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1. Introduction 

Impurity ions present in tokamak plasmas can significantly degrade plasma performance 
through radiative losses and fuel dilution [1]. Conversely, impurity ions present in the 
edge or near the divertor of a tokamak may be beneficial by radiating thermal energy and 
thereby mitigating the heat load incident on the divertor and plasma facing components 
[2]. An understanding of impurity transport (a topic which has received less attention 
than the transport of bulk ions) is crucial in order to predict the overall performance of 
fusion plasmas [1]. 

The analysis of charged particle motion in strongly-magnetized plasmas, for 
example those in tokamaks, generally relies heavily on a multi-timescale approach [3], 
whereby the particle motion is decomposed into fast gyromotion about the magnetic 
field, streaming motion along the field (characterised by adiabatic invariance of the 
particle magnetic moment) and a much slower drift motion of the particle guiding centre 
across the field. This treatment is valid provided that the drift ordering is satisfied. 
Specifically, the particle Larmor radius must be small compared to the length scales on 
which the electric, E, and magnetic, B, fields vary, the cyclotron period must be short 
compared to the field variation time scales and the E x B drift velocity must be small 
compared to the particle thermal speed. 

In spherical tokamaks (STs), such as the Mega Ampere Spherical Tokamak (MAST) 
[1] , the drift ordering is often not applicable, and it is then not possible to treat particle 
motion analytically. On the low field side of a ST plasma the poloidal component of 
the magnetic field can be comparable to the toroidal component, so that the drift orbit 
width of a trapped particle can be of the same order as the particle Larmor radius. 
Moreover, in ST plasmas with edge transport barriers (ETBs) the associated radial 
electric field can vary on length scales comparable to both the drift orbit width and the 
Larmor radius of thermal ions 0]. 

It has been shown that both coUisionless trapped particle guiding centre orbits and 
particle gyro-motion are significantly distorted by strongly-sheared radial electric fields 
O E]. This effect is often referred to as orbit squeezing [7], although in fact the widths 
of trapped particle orbits can be either reduced or increased, depending on the sign 
of the electric field gradient. The effects of orbit squeezing on neoclassical bulk ion 
transport have previously been studied in the limit of large aspect ratio [5], |8] using this 
framework. However, the analytical theory of orbit squeezing assumes linearly sheared 
electric fields, and recent measurements of radial electric fields associated with ETBs in 
spherical tokamaks indicate field profiles that are not consistent with the assumption of 
constant shear |4j. This highlights one limitation of the current theory and motivates 
the use of a numerical approach. 

We note that non-uniformly sheared electric fields have previously been investigated 
using a full orbit approach [9]. However, this study investigated the impact of shear 
on particle E x B flow and only considered coUisionless particle orbits. In the present 
work we investigate the impact on the motion and collisional transport of impurities of 
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a strongly-sheared electric field in a MAST-like spherical tokamak equilibrium using a 
full orbit, test-particle simulation code. The fields are assumed to be axisymmetric and 
static; we therefore neglect turbulent transport effects. 

2. Model 

The transport of impurities in a spherical tokamak with a prescribed local radial electric 
field is simulated using the test particle code CUEBIT [IHIIII], which solves the Lorentz- 
Langevin equation 

dv m 

m^— = Ze (E + V X B) ^ (v — u) + m^a^, (1) 

dt r^i 

where m^, Ze and v are the impurity mass, charge and velocity, respectively. The 

term on the right-hand side of ([T]) that is proportional to v — u models the friction 

from the bulk ions. In the absence of an electric field it ensures that the impurities 

relax to a drifting Maxwellian distribution temperature Ti and flow u matching those 

of the bulk ions. Here, u is taken to be zero; we assume that the radial Lorentz 

force is sub-dominant to the pressure gradient term in the radial component of the 

bulk ion fluid momentum balance equation. This assumption is consistent with 

reported measurements of temperatures, densities and flows in the vicinity of ETBs 

in MAST [U [12]: the pressure gradient is typically several tens of kPam~^ whereas the 

contributions of poloidal and toroidal flows to the radial component of the Lorentz force 

are typically around one or two kPam~^ at most. In taking u to be zero we also neglect 

any deviation of the bulk ion distribution from a stationary Maxwellian arising from 

neoclassical effects [13]. The quantity r^i is the classical impurity-bulk ion collision time 

given by [13] 

which is a function of the local bulk ion temperature and density profiles, T, and rzj, 
respectively. Here, In A is the Coulomb logarithm and eo is the permittivity of free space. 

The final term on the right-hand side of ([1]) models random accelerations of 
impurities due to Coulomb collisions with bulk ions. The vector is a set of three 
random numbers, chosen independently for each particle at each time step from a 
Gaussian distribution of zero mean and variance 

^ = ^t^ 

where At is the time step that is used in the numerical simulation. Although the 
instantaneous collisions are taken to be isotropic, the cumulative effect of these over 
an impurity Larmor orbit naturally reflects the effects of gradients in the bulk ion 
temperature and density profiles (leading to, for example, thermal forces), which leads 
to anisotropy in the collisions. Here we have again assumed that the bulk ions have a 
Maxwellian distribution. 
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Furthermore, in simulating test particles we assume that Z'^n^ <^ ni (me/mi)^^^, 
where Z is the impurity ion charge number, is the impurity ion number density, rij 
is the bulk ion number density and me and mj are the electron and bulk ion masses 
respectively [13]. This condition essentially ensures that collisions of bulk ions with 
impurity ions are sufficiently infrequent that they have a negligible effect on the bulk 
plasma neoclassical transport and automatically ensures that collisions between impurity 
ions can be neglected. 

A MAST-like equilibrium is considered, which is modelled using an analytical 
Solov'ev-type solution of the Grad-Shafranov equation with a linear pressure profile 
and a potential toroidal magnetic field component. In cylindrical coordinates (R,0,Z), 
the poloidal flux, is of the form 

+ ^^R'z'Y (4) 

with Rq = 0.964 m, i?b = 0.93 m, 7 = 0.8 and \E'o — 0.9 Tm~^ corresponding to 
a total plasma current of 1 MA. The magnetic field vector is calculated from (jl]) as 
B=i?i?<^V0+V^' X V0 where RB^= 0.386 Tm. We assume Ti and rii are flux functions 
and prescribe analytical forms, 

Ti = To-5^ + Ti, (5) 

ni = no— + ni, (6) 

where \E'm the flux at the magnetic axis and Tq = 1 keV, Ti = 0.5 keV, no = 5 x 10^^ 
m~^ and rii = 1 x 10^^ m~^. These profiles are broadly consistent with those measured 
in co-rotating MAST H-mode plasma discharges In such discharges, the effective 
ion charge is often measured to be close to unity throughout the core plasma, indicating 
no significant impurity accumulation. 

A local radially sheared electric field E = — V$ near the edge is included, with an 
associated electrostatic potential $ of the form 

/\E' — \E'i \ 

$ (^) = $0 arctan — — — ^ , (7) 



where $0? ^1 and A\Ef are the potential barrier height, location in space and width 
in space, respectively. Measurements reported by Meyer et al [4] of edge radial 
electric fields in MAST H-mode plasmas (see for example, figure 7 of ^) show that 
the corresponding potential proffie is approximately of the form given by ([7]) and that 
these constants have typical values of $o=60 V and A\E'=1.35 10"^ Tm^. We choose, 
for computational ease, to centre the electric field structure on \I'i=0.6\E'(-Ro, 0)=- 
0.4 Tm^. For our particular choice of equilibrium, these quantities correspond to a 
physical location and width in the outboard Z=0 plane of i?i=1.2 m and Ai?=l.l 
cm. The corresponding peak electric field strength, Eq, is —10.5 kVm~^, which is close 
to measured values in MAST ETBs We shall present results from simulations 
with varying values of $0 a^^d A\Ef within the ranges —4.75 < $0 ^ —120 and 
0.01\E'(i?o,0) < A\E' < 0.08\I'(i?o,0). The form of this prescribed electric field is shown 



Collisional impurity transport and sheared electric fields 



5 



in figure [H which illustrates a limitation of using the orbit squeezing model: namely, 
which value of the squeeze factor to use when characterising the impact of the electric 
field on particle motion. 

In addition to a radial electric field, the inward Ware pinch [15] is modelled 
approximately with a constant toroidal electric field = —0.3 Vm~^. Including 
the 1/R dependence of this field has a small, 0(< 10%), effect on measured diffusion 
coefficients. 



3. Characterising transport 

In the absence of an electric field we expect neoclassical particle transport to be purely 
diffusive [13]. For tokamak plasmas with slowly- varying profiles, transport coefficients 
may be deduced empirically from moments of the spatial distribution of test particles 
[T6| [T7] . Figure [21 in which the evolution of the mean minor radius of the impurity 
distribution is plotted, illustrates a typical CUEBIT simulation in which 10^ C^"*" ions are 
initially released from the magnetic axis in the absence of any electric field. Overlaid is a 
t^/^ fit expected from diffusive transport [18]. This corresponds to a diffusion coefficient 
of i5 ~ 10 m^s~^, which is close to the expected Pfirsch-Schliiter neoclassical diffusivity 
in the plasma core, D ~ (fpt/Tzi- Here q is the safety factor and pl is the Larmor 
radius. 

In the case of narrow transport barriers, with plasma properties varying on length 
scales down to the Larmor radius, it is not possible to infer local transport coefficients 
using this method. In such local effective diffusivity D^q may be deduced directly 

from the local radial particle flux and density gradient dn^/dp: 

dn^/dp' 

Here is the flux surface- averaged minority ion density and p = p(^I/) is a radial 
coordinate, defined as the flux surface averaged minor radial coordinate 

= / /' d<i/', (9) 



/ Rdl 

with / the arc length along a flux surface in the R-Z plane. p(^') is evaluated numerically 
as a function of for the Solov'ev equilibrium described above. To quantify the impact 
of sheared electric fields on particle transport the value of D^s is normalised to its value 
measured in the absence of any electric field, D. 

Clearly is a simplified formulation of the particle diffusion coefficient since we 
assume that particle flux scales linearly with the minority ion density gradient. However, 
we wish to emphasise that Fz, which appears in ([S]), is the full particle flux as computed 
by CUEBIT. This includes the effects of transport driven by bulk ion gradients, e.g. the 
thermal force arising from the bulk ion temperature gradient (see, for example, equation 
5.9 of [13]). Such a force arises from the variation of the collision frequency across a 
Larmor orbit and, unlike to guiding centre and fluid calculations, appears naturally in 
full orbit particle simulations. 
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In order to elucidate the impact of radially sheared electric fields on particle 
transport we use ([1]). In order to proceed analytically we neglect both the random 
Langevin term and the bulk ion fiow velocity, whose contribution to the momentum 
balance equation can, as noted previously, be assumed to be sub-dominant to the bulk 
ion pressure gradient term. The combined effect of the electric field and collisional drag 
terms can be seen by adopting Cartesian coordinates, with B = i?^ and E = both 
assumed to be constant, where we use the x-direction as a proxy for the radial direction 
in tokamak geometry. Equation ([T]) can then be solved exactly, yielding 

v^ = v±exp {-ly^it) sin llt + — — ^-^-^j^, (10) 

Vy = v±exp{-u^it)cosQt - ^Y^-j^y^. (11) 

In tokamaks we have that z/^i ^ f^^. In this limit, we are left with the usual E x B drift 
in the y-direction and an additional "drag force" drift in the x-direction, which arises 
from the inclusion of a collisional drag in the test particle equation of motion. This 
implies an inward-pinch in the presence of a negative radial electric field. Substituting 
for the impurity- ion collision frequency, given by ([2]), we find that the drag force drift 
can be written as 

" ^F5f75^^|-^- (12) 

This is independent of test particle mass but proportional to test particle charge. 

In order to relate the impact of a collisional drag force drift on particle transport 
we consider the steady-state one-dimensional transport equation 

_d_ ijj^ _^ 

dx I dx 1 + (x/Ax)2 

where v is the peak value of the drag drift, v^- Here we have assumed a form of the 
spatial variation of E^ that is consistent with the electrostatic potential given in (^^, 
with X the displacement from the peak electric field (x = ^ — ^i) and Ax a measure 
of the electric barrier width (Ax = A^^). We neglect all spatial variations in Vd except 
for that occurring due to its dependence on E^. We also assume that the impurity ion 
Larmor radii are small compared to Ax and that particle drifts in the x-direction are 
due solely to the drag effect. The above form of the transport equation ensures that 
the test particle flux, which we assume is composed of the usual diffusive and advective 
terms, 

T, = -D^-Vdn,, (14) 
dx 

is constant. Introducing a dimensionless spatial variable ^ = x/Ax and a Peclet number 
Pe = vAx/D, ( TT^ becomes 

dn^ Pen^ 



0, (13) 
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where 7 = r^Ax/D. This has the exact solution 



r 

(0 = "-0 exp Petan^^ 1 — j exp (^Pe tan~^ r/^dr/ 



(16) 



where uq is the particle density at a; = 0, i.e. at the position of peak electric field and 
drag drift. For ^ ^ 1, or x ^ Ax, this reduces to 



n^{x) = noexp (^-^^) - 
It follows that, at X = 0, we have 

r. D 



V 



exp 



Pex 



Ax / 



(17) 



driz/dx 1 + vno/Tz 

Identifying v with the drag force drift velocity, and assuming that n^/T^ is essentially 
constant between simulations, we therefore find that the local effective particle diffusion 
coefficient scales with test particle charge as 

D 



D. 



eff 



where 



a 



1 + aZ' 
Ex n^Ui m}J'^e^ In A 



(19) 



(20) 



and is independent of both the test particle mass and charge. Despite the 
approximations used to derive (1191) . we will show in the next section that it provides a 
fairly accurate description of impurity ion transport at an ETB in full toroidal geometry. 



4. Impurity transport simulations 

4.1. Transport scaling with particle parameters 

We now present simulation results for a number of impurity species. In table [H we list, 
for reference, the normalised collision frequency. 



of each species at the position of the simulated transport barrier. Here f th is the particle 
thermal velocity. The quantity v* measures the number of collisions a particle undergoes 
in one toroidal transit of the tokamak, with u* < 1 indicating that a species is in the 
banana regime and u* > 1 that a species is in the Pfirsch-Schliiter regime. 

Figure |3] shows the results from simulations in which 10^ ions of various impurity 
species were released from the magnetic axis in equilibria with and without an inward- 
directed, radially sheared electric field. The electric field acts to confine C^"*" ions. The 
comparison of density profiles of C^"*" [see dashed lines, figure |3l^a) and figure |3]^b)] shows 
a step in the density profile at a minor radius of approximately 0.4 m. This indicates 
that the sheared electric field is acting as a barrier to the outward radial transport of 
the ions. The location of this step coincides with the peak of the radial electric field. 
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Table 1. Normalised collision frequency for impurity ions simulated with computed 
effective diffusion coefficients, ([5]), without, Dq, and with, DcS, a sheared radial electric 
field as given by ([7]). For those simulations indicated by * we have modified the location 
of the barrier, such that ^I^i = 0.3^'(i?o,0) — —0.2 Tm^, and temperature and density 
profiles, such that Tq = 0.75 keV, Ti = 0.05 keV, riQ = 5 x 10^^ m^^ and rii = 2 x f 0^^ 
m . 



Impurity species 


l/* 


^0 




He2+ 


0'^2 


29.7 ± 1.4 


± ± . t _i_ u .0 




0.f67 


7.8 ±0.5 


0.8 ±0.2 


Ne2+ 


0.0f4 


43.7 ±0.8 


33.0 ±0.5 


Ne4+ 


0.057 


18.5 ±0.6 


10.0 ±0.3 


Ne^+ 


0.129 


10.7 ±0.7 


4.1 ±0.2 


Ne«+ 


0.229 


6.3 ±0.4 


2.2 ±0.2 


Neio+ 


0.358 


4.5 ±0.3 


1.5 ±0.1 




4.737 


2.8 ±0.3 


0.75 ±0.07* 


Ario+ 


0.253 


5.3 ±0.3 


1.6 ±0.1 




3.349 


2.7 ±0.2 


0.49 ±0.05* 


M0IO+ 


0.167 


6.4 ±0.8 


1.7 ±0.08 




2.208 


3.4 ±0.2 


0.36 ±0.03* 




0.118 


6.2 ±0.2 


2.2 ±0.07 




1.561 


3.5 ±0.1 


0.29 ±0.02* 


•^•20+ 


0.473 


1.7±0.1 


< 0.2 



The variation of impurity density with minor radial coordinate for three other 
impurity species, He^"*", Ne^°"^ and W^""*", are shown on figure O The density profiles 
are similar in the absence of an electric field but differ substantially when it is present. 
The greatest confinement occurs for W^°"^. All impurity species simulated have reached 
thermal equilibrium with the bulk ions particles before they begin interacting with the 
radial electric field. The poloidal distribution after 100 ms of simulated time are shown 
in figure m for He^"*", Ne^""*" and W^""*" ions, which again clearly illustrates the enhanced 
confinement of impurity ions by a sheared radial electric field, particularly in the case 

of W20+. 

In order to compare the scaling of effective diffusivity with particle parameters, 
given by f|T9|) . with simulations we consider results for a number of test particle species 
with varying mass and charge numbers. In each case the simulation is run until the 
ratio given by ([8]) at the peak of the electric field has reached a steady-state, except for 
noise fluctuations arising from the use of a finite number of test particles. However, the 
absence of a continuous particle source means that both the particle flux and the density 
gradient continue to slowly decline. In figure [5^ the normalised diffusion coefficient is 
plotted for particle mass numbers in the range 20 (neon) to 184 (tungsten) all with 
Z=10. We see that the diffusion coefficient is essentially independent of mass number. 
In figure [5b the test particle mass number is held fixed at 20 but the charge state is 
varied from Z=2 to Z=10. Equation (JTOil predicts that the quantity D/Des — 1 varies 
linearly with particle charge number. This is broadly consistent with the results shown 
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in figure 03, although the uncertainties in D/D^s — 1 increase with Z. Furthermore, in 
deriving f|T9|) contributions to the particle flux from bulk ion density and temperature 
gradients were neglected (cf. f |T^ and equation 5.9 of pjj). Including such terms would 
add to (fT9l) a term that is non-linear in Z in its denominator and which could explain 
the possible turn-over present for high Z in flgure [5l 



4-2. Transport scaling with electric field strength and width 

Figure E] shows the variation with electric fleld width and height of the effective diffusion 
coefficient for C^"*" ions. The normalised effective diffusivity has a power law scaling with 
electric barrier width for constant barrier height. The coefficients of the fit themselves 
scale linearly with barrier height. Thus the scaling of diffusion coefficient with electric 
fleld width and strength may be expressed in the form 

where AR is the FWHM of the electric fleld in the outer midplane and E is the peak 
height of the radial electric fleld barrier. We see that the diffusion coefficient scales 
as the width of the electric fleld barrier to a power of the electric fleld strength, with 
the width of the barrier normalised to a fraction of the bulk ion Larmor radius. This 
scaling has been verifled for AR and E in the ranges 0.5pu < AR < Apu and —0.65 
kVm~^ < E < —7.65 kVm~^ for all impurity species, with pu the bulk ion Larmor 
radius. In the limit of no electric fleld, this scaling recovers the effective diffusivity as 
found in simulations without an electric fleld barrier. The scaling parameters for various 
species are listed in table El 

Unlike (fT9|) . the empirical expression for D^s given by (122|) does in fact depend on 
the barrier width, implying limits in the validity of the approximations used to derive 
the former. We expect that for large AR/pu, the effective diffusion coefficient will no 
longer have a power law scaling but will tend asymptotically to a result similar to ( TT9l) . 

We note, from table [21 that ZEq and ARq are broadly constant. Therefore, the 
dependency of impurity transport on the shape of the electric fleld barrier depends 
approximately only on the charge state of the species, not on the species mass. Thus, 
we can approximate this scaling as 

Deff ^ Do — ^ . (23) 



0.2p 



Li , 



5. Discussion 



We have simulated realistic radial electric fleld proflles for ETBs in a spherical tokamak 
equilibrium and investigated the impact on collisional impurity transport of the strength 
and width of the electric fleld proflle for a variety of impurity species. We have 
demonstrated that full orbit simulations provide a valuable insight into the transport 
of particles in such scenarios. It has been shown that a strongly-sheared radial electric 
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Table 2. Scaling of computed effective diffusion coefficient, (jS]), with electric field 
strength, E, and width, Ai?, for several impurity species. 



Impurity species 


Do Ks-i) 


Ai?o 






Eo (kVm-i) 




8.3 ±0.4 


(0.17 


± 


0.02) pLi 


8.1 ±0.4 


Neio+ 


3.4 ±0.4 


(0.21 


± 


0.03) pLi 


5.7±0.4 




2.0 ±0.6 


(0.20 


± 


0.02) pLi 


4.2 ±0.4 


Mo2"+ 


f.6±0.8 


(0.22 


± 


0.03) pLi 


2.9 ±0.6 



field can significantly increase the confinement of test particle ions. This effect has been 
explained in terms of a drift arising from the drag between impurity and bulk ions, with 
the direction of the drift given by the radial electric field and its strength determined 
by the impurity charge. Clearly, the negative implication of this study is the enhanced 
confinement of highly ionized impurity species in the vicinity of strong, negative radial 
electric fields. 

In our derivation of the drag force drift given by (fT2|) we neglected both spatial 
variations in the radial electric field and toroidal geometry. However, the scaling of 
effective diffusivity with Z inferred from using this drift velocity in a simple slab model 
of the transport barrier, ( TT9l) . is in agreement with our test-particle simulations in which 
no approximations are made with regard to finite Larmor radius or toroidal effects. 
This suggests that test-particle transport in the immediate vicinity of the barrier is not 
strongly-dependent on the effects of toroidal geometry or the collisionality regime (i.e. 
the value of z/*). This further indicates that finite Larmor radius effects are relatively 
unimportant. One would expect electric field shear on length scales approaching the 
particle Larmor radius to result in a modified drag force drift velocity, analogous to the 
E X B drift with finite Larmor radius correction discussed by Tao et al [9]. However, 
this requires further study. 

The location of the electric field in the simulations, which resembles more an internal 
transport barrier than an ETB, was chosen in order to lower computation time by 
reducing the distance from the plasma core to the location of the electric field, thereby 
reducing the time before particles begin interacting with the electric field. This also 
has the benefit of improving particle statistics, since the particles are less dispersed 
radially and the number of particles at the barrier is consequently higher. Simulations 
in which the electric barrier is located closer to the plasma edge, listed in table [H 
demonstrate that moving the barrier closer to the edge has little qualitative impact on 
the results obtained. Furthermore, the variation of collisionality of impurity species, 
with normalised collisionality ranging from 0.014 (banana regime) to 4.74 (Pfirsch- 
Schliiter regime), appears to be unimportant, since the measured normalised diffusivity 
remains independent of mass and scales as expected from the drag force drift argument 
with particle charge. 

Whilst the effect of bulk plasma rotation was neglected in the present study, one 
would expect transport barriers to be associated with flows of the bulk ions (e.g. zonal 
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flows). For impurity ions in particular inertial effects will become important in rotating 
plasmas. As recently shown [17] , strong toroidal rotation can cause heavy, incompletely- 
ionised impurities such as W^°+ to undergo very rapid collisional transport (due to the 
combined effect of trapping in a centrifugal potential well and a modification to the 
effective magnetic field arising from the Coriolis force). One would therefore expect 
a competition between this enhanced transport due to rotation and the improved 
confinement resulting from sheared radial electric fields. 

We note that the effect of improved confinement due to a sheared radial electric 
field is more important for spherical tokamak plasmas than for conventional tokamaks, 
primarily due to the strong B dependence in ( JT2l) . For example, in the MAST-like 
equilibrium used throughout, we have that the Ware pinch velocity is of the order 
vware — qE^/eB ~ 5.6 ms~^ and the drag force drift velocity is of the order Vd ~ 8.2 
ms~^ for C^'^ (where we have used B = 0.4 T, Ex = —10 kVm~^, E^ = —0.3 Vm~^, 
rzj = 3 X 10^^ m~'^, Tj = 500 eV, q = 3 and e = 0.4). Typical values for a large 
conventional tokamak such as JET, on the other hand, are fware ~ 0.4 ms~^ and 
Vd ~ 0.02 ms-^ (assuming 5 = 3T, = 20kVm-\ E^ = -0.04 Vm-\ = 3 x 10^^ 
m~^, Tj = 5 keV, q = 3 and e = 0.1 [19]). Thus, the drag force drift is much less 
significant than the Ware pinch under typical conventional tokamak conditions. 

An area of possible future investigation is the impact of fiuctuating radial electric 
fields on collisional particle transport, with the fiuctuations resulting from, for example, 
geodesic acoustic mode (GAM) oscillations [20]. Furthermore, in recent fiuid simulations 
of the COMPASS-D tokamak, it has been observed that the H-mode radial electric 
field close to the plasma edge changes from being negative between edge localised 
modes (ELMs) to positive during an ELM [21]. In the context of the present work 
one would expect the direction of the drag force drift to switch between inward and 
outward directed, which could have important implications for impurity accumulation 
and the transport of particles at transport barriers. The fiexibility of the full orbit 
approach means that turbulent electric or magnetic fields could be readily incorporated 
into numerical studies of impurity transport in ELMy H-mode plasmas. 
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Figure 1. Form of the prescribed electric field (solid line) and corresponding squeeze 
factor for a C^"*" (dashed line) ion in the outer midplane of a MAST-like equilibrium. 
The full-width half-maximum of the radial electric field is chosen to be equal to the 
bulk ion gyroradius evaluated at the location of maximum electric field. The squeeze 
factor is calculated from 5 = 1-1- Ze'^"{r)/m^Vl^. 
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Figure 2. Mean minor radial coordinate versus time of C^''" ions released at the 
magnetic axis in a simulation with zero electric field (squares). Overlaid is a fitted i^/^ 
curve expected of diffusive transport (solid line) . 
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Figure 3. Density profiles of impurity ions against minor radius, both without (a) and 

with (b) a sheared radial electric field. Simulated ions are He^+ ( ), C^"*" ( ), 

]\[gio+ ^ — . — ^ g^jjj^ ^20+ ^ y Vertical lines in (b) indicate the peak strength and 

full- width half-maximum positions of the sheared radial electric field. 
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Figure 4. (Colour online) Poloidal distribution of (a) He^+ (b) Nei°+ and (c) W20+ 
ions in a MAST-like equilibrium with a radially sheared electric field, the position of 
which is indicated by a dashed (red online) line. 
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Figure 5. Scaling of particle diffusion coefficient with particles parameters, (a) scaling 

with test particle mass with constant charge number Z=10. (b) scaling with test 
particle charge for constant mass number of 20 (neon). 
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Figure 6. Contours of constant normalised effective diffusion coefficient for ions 
with radial electric field height and width. 



